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Abstract 

We present evidence for the existence of infinitely-many new families of renor- 
malisation group flows between the nonunitary minimal models of conformal field 
theory. These are associated with perturbations by the 021 an d (f>i5 operators, and 
generalise a family of flows discovered by Martins. In all of the new flows, the finite- 
volume effective central charge is a non-monotonic function of the system size. The 
evolution of this effective central charge is studied by means of a nonlinear inte- 
gral equation, a massless variant of an equation recently found to describe certain 
massive perturbations of these same models. We also observe that a similar non- 
monotonicity arises in the more familiar <fii3 perturbations, when the flows induced 
are between nonunitary minimal models. 
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1 Introduction 



Integrable quantum field theories which interpolate between different conformal field 
theories have been recognised as an intriguing feature of two-dimensional models ever 
since the initial paper of A.B. Zamolodchikov O. The first examples arose as 0i3 
perturbations of unitary minimal models M PjP+ i. The existence of higher-spin integrals 
of motion Q indicates that such perturbations should be integrable; in Q (see also j|) 
perturbative arguments were used to show that for one sign of the coupling constant 
the resulting flow is to the neighbouring unitary minimal model, namely M p -i tP , at 
least for p > 1, Further support for this picture, this time valid at all values of p, 
was provided when ideas from the thermodynamic Bethe ansatz (TBA) enabled exact 
integral equations describing the flows to be proposed p). 

Unitarity of a conformal field theory is not a necessary requirement for the existence 
of integrable perturbations, and in |||6| the arguments of |l| were extended to the 0i3 
perturbations of more general minimal models Mp, q - The results indicated that there 
should exist integrable flows between the models A4 Ptq and M.2p-q,p, at least in the 
region p 3> q—p where calculations perturbative in {q—p)/p can be trusted. (It remains 
an open problem to give a TBA treatment of the nonunitary flows, which would not 
suffer from this caveat.) The unitary and nonunitary models can be put on a common 
footing by defining a parameter ( = p/ (q— p)Q The value of £ suffices to identify p and 
q uniquely, since they must be coprime. In all cases the predicted flow is from the model 
specified by ( to the model specified by (— 1. 

Many further massless flows have since been discovered and studied, often by means 
of the TBA technique: the papers f7Hl9| are a sample of this work. Our main interest 
here will be the flows found by Martins and further studied by Ravanini et al. fl9[| . 
These are perturbations of the minimal models M. v ^,p-\ and M. Pj 2p+i by the operators 
02i and respectively, and the predicted trajectories are: 

02i : M Pt 2 P -i -> M p -i,2p-i ; (1-1) 

015 : M P) 2p+l -> M Pt 2p-l ■ (1.2) 

These flows can be chained together to form a single sequence, along which the perturb- 
ing operator alternates between 02 1 and 0i5. 

In this paper we will extend this picture by proposing further sequences of flows for 
which the perturbing operator alternates in the same way. As in studies based on the 
TBA technique, our main tool will be a conjecture for an exact equation expressing the 
finite-volume ground state energy of each model as a function of the system size R. In 
fact, it will be most convenient to work with an associated scaling function known as 
the 'effective central charge', c e ff(r). This is related to the ground state energy as 

E (M, R) = E hulk (M, R) - ^c cfr (r) , r = MR . (1.3) 

was denoted m in 0; the symbol p is also often used, but we are reserving this for the integer- 
valued first index of the minimal models M P , q 
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In this equation, Eq is the ground state energy, -B^uik the irregular bulk part, and M some 
mass scale, set for example by an infinite- volume one-particle state. The effective central 
charge depends on the system size only through the dimensionless combination MR; for 
a model which happens to be scale invariant (conformal), it is constant and equal to 
c— 24cZ, where c is the model's central charge and d its lowest conformal dimension. 
The previously-studied sequence ( |l . l|) , (|l.2| ) is picked out by the monotonicity of the 
function function of r. In all other cases, somewhat to our surprise, we found 

that the effective central charge undergoes a number of oscillations as it interpolates 
between its short and long distance limits. 

Our approach differs from the TBA method, and is much closer to that of the pa- 
pers [p0|-p6|, in that a single nonlinear integral equation (NLIE) is proposed to describe 
infinitely-many different perturbed conformal field theories, each being picked out by an 
appropriate choice of certain parameters. For massless flows, the one previous example 
of such an equation was found by Al. Zamolodchikov |^] , and his paper formed a large 
part of the motivation for our work. 

The particular nonlinear integral equation that we will be using is obtained in §2; 
then in §3 we discuss the nature of the flows that it predicts, and in §4 we take a 
more detailed look at various asymptotics. This allows us to back up our conjectures 
with a comparison with UV conformal perturbation theory. In §5 we comment on some 
features of the </>i3 flows, and in the concluding §6 we indicate some open problems that 
remain for future work. 



2 The nonlinear integral equation 

The work of [20] built in part on the nonlinear integral equation of [21]], which encodes 



the finite-volume ground state energy of the sine-Gordon model for a continuous range 
of the coupling (3, and at general 'twist' a. Setting ( = [3 2 / (8-k — (3 2 ) , it can be shown 
that for ( = p/ (q—p) and a = 1/p, the ground-state energy of the 0i3 perturbation of a 



general minimal model Mp,q is also matched [23,24, 27|-|29|] . (Note that this definition 
of Q is thus in line with the one given in the introduction.) Since the sine-Gordon 
model is massive, the minimal model perturbations reproduced in this way are also 
massive. In contrast, the modification to the equation found in |2(| describes a massless 
perturbation of the sine-Gordon model, flowing from the model at £ to the model at 
£— 1, and thus matching the pattern of massless <j)\z flows discussed in the introduction. 
While only the a = case was treated explicitly in |2(|, we will confirm below that with 
a suitable nonzero value of a the massless ^13 perturbations of minimal models are also 
obtained, with results that agree with the perturbative picture. In the nonunitary cases 
this yields some genuinely new information, since, as already mentioned, TBA systems 
giving the exact evolution of the ground-state energy for these models are not known. 
However this is not the main purpose of our paper: rather, we would like to perform a 



Hhe scale-invariant form of this equation had arisen previously in the context of integrable lattice 
models El 
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similar trick for the 4>2i and 0is perturbations of minimal models. Recently, a nonlinear 
integral equation describing finite size effects in the a 2 model was conjectured [25 1], 



and checked to describe the massive c/>i2 , ^21 and </>i5 perturbations of minimal models 

(2) 

on the imposition of a suitably-chosen twist. (The connection between the a 2 model 



and these perturbations can be understood through quantum group reduction |30|-|32j.) 
Our aim is to modify this equation so as to find out about the massless, interpolating, 
flows that may also be induced by these perturbations. 
The nonlinear integral equation found in |2f| is 

f(p) = ina- ir sinhfl + / tp(0-e') ln(l+e /(e,) ) d6' - [ <p(6-0') ln(l+ e - /(e,) ) d6' (2.1) 

Jd Jc 2 

with the effective central charge given in terms of its solution as 



c e s(r) 



^j(^J smh8ln{l+e f{e) )d6 - J sinhfl ln(l+ e - /(9) ) d6^j . (2.2) 



The contours C\ and C2 run from —00 to +00, just below and just above the real #-axis, 
and r is equal to MR, R being the size of the system and M the mass scale, set by the 
fundamental kink. The kernel 

y°° e ifcfl S mh(p)cosh(gfc(l-20) dfc 
n ) cosh(ffc)sinh(A;fC) 2vr ( ' ' 

is equal to i /2n times the logarithmic derivative of the scalar factor in the S-matrix of the 

(2) 

a 2 (Tzitzeica-Izergin-Korepin-Bullough-Dodd-Zhiber-Mikhailov-Shabat . . . ) model. 

(2) 

The parameter £, related to the a 2 coupling 7 as £ = j/(2tt — 7), is the analogue 
of the sine-Gordon parameter £ mentioned earlier, while a corresponds to the twist. To 
obtain massive 0i2, 4>2i an d 1^15 perturbations of a minimal model -A4p j(? , the values of 



£ and a must be chosen as follows f25fl : 



>12 



'>21 



i = ^— , a = 2/p , p<2q; (2.4) 







1 


p 






1 




■2p 




1 


<-l 






1 





i = 2 —— , a = 2/q , p >q/2 ; (2.5) 



2p 1 



In each case, the inequality delimits the region of the (p, q) plane within which £ is 
positive, and the corresponding perturbation is relevant. (The first inequality is not 
strictly necessary, since we will be adopting the convention that p < q in the specification 

*as for the sine-Gordon model, a scale invariant version of this equation had previously been found 
in the context of integrable lattice models [M 
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of A4 P; q.) The UV effective central charge, c e g(0) = 1 — ^pj-a 2 , is always equal to 



1 — 6/pq, but the subsequent terms in the expansions differ, as expected given the 
different perturbations being described. 



The modification to the massive equation of [21] found in [2C] made use of elements 



of the corresponding massless scattering theory, proposed in [33]. In our case we do not 
have a description of the massless scattering in terms of an S-matrix but, proceeding 
by analogy with the results of pcfl , we will substitute the single equation ( |2.1|) with 
two equations describing hypothetical left and right movers. Prom the formulae ( J2.5D 
and ( |2.6[) , we notice that the massive versions of the perturbations (|1 . 1| ) and (|1.2| ) have 
£ = 2p — 1 and £ = 2p respectively. We therefore seek a massless modification of 
the NLIE ( |2.1| ) which will interpolate in general between an ultraviolet theory with 
parameter £ and an infrared theory with parameter £ — 1. To this end, we introduce 
two analytic functions /r(0) and fh{9), couple them together via 



— i—e +VKa 

+ [ 0(6 -e')ln(l + e fR(e ">) dO'- [ <f)(9-9')ki(l + e- fR ^) dd' 
Jd Jc 2 

+ [ x{0-0')\n{l + e-f L ^) dd' - f x(0-0')Hl + e fL{e,) ) d9' (2.7) 



hie) 



-i— e — nra 



+ / 0(0-0') Ml + dO' - / 4>{6-6')\n{l + e- fL{el) ) d9' 

IC 2 JC! 



+ [ x{0-0')\n(l + e- fR ^) d9'- [ X {0 ~ 0') ln(l + e fR ^) d9' (2.8) 
Jc 2 JCx 

and replace the expression (|2.2|) for the effective central charge with 



3ir 

Ccff(r) = 2^ 



/ e e ln(l+e /fl(e) ) dO - [ e e ln(l+ e - /fl(9) ) dO 
JCi Jc 2 



+ 



'ln(l+e 



d£ 



Cz 



'ln(l+e 



d,0 



Ci 



As before, the contours C\ and C2 run from — 00 to +00 just above and just below the 
real axis, and r = MR. However, since the flows are massless M no longer has a direct 
interpretation as the mass of an asymptotic state, but rather sets the crossover scale. 

In the far infrared, r — > 00 and the two equations ( |2.7D and ( |2.8| ) decouple. The 
result is two copies of the UV limit of (|2.l|) , with the kernel tp{9) substituted by 4>(9). 
Since the IR destination is to be the model at £ — 1, we take 4>(9) to be the massive 
kernel ( |2.3[) with the substitution ^ — ^ ^ — 1: 



00 gik9 



smh(fcf ) eosh(|A;(3-2£)) dk 



cosh(f A;)sinh(/cf (£-1)) 2vr 



(2.10) 
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To find the other kernel, x(#), we consider the UV behaviour of the new system. 
This should coincide with the UV limit of the massive system (|2.1| ). As r —* 0, the 
solution to fl2.7P , ( |2.8| ) splits in the standard way into a pair of kink systems, centered 
at 9 = ±ln(l/r). We focus on one of them by replacing 9 by 9 — ln(l/r), and keeping 
this variable finite as the limit is taken. The exponential term in ( |2.7| ) is replaced by 
— T^e® , while to leading order the exponential in (|2.8| ) can be neglected. Furthermore, 
the singularities in ln(l + e^^ L ^') previously found on the real axis are pushed to — oo, 
allowing integration contours in any integral involving fa to be moved across the real 
9 axis at will. It will be convenient to shift C2 down to coincide with C\ in all such 
integrals, and to take C\ to be the line 9m 9 = —77, with rj a suitably-small real number. 

For /r(9), the singularities do not disappear but using the reality properties of Jr we 
can at least rewrite the integrations above and below the axis in terms of the imaginary 
part of a single integration along C\. The final equation is 



00 



f R {9) = -!-e d + ma' + 2i / 0(0 - (9' - ij7))9fm(ln(l + e^"^)) d9' 



3C 



x(e-(9> -iv))fa(9' -in) d0' (2.11) 
fa{9) = -W + / ^{9-{9' + i v ))f L {9' + d9' 

/oo 
X {9 - (9' - M7))3m(hi(l + e fR(0 '- iv) )) d0' . (2.12) 
-00 

Taking the Fourier transform of ( |2.12| ) yields 

(1 - 4>(k)) fa(k) = -2i 7r 2 a'S(k) - 2 % x(k) L R (k) (2.13) 

where Ln(k) denotes the transform of Qm(ln(l + &f R ^~ ir ")) and our convention for the 
Fourier transform of f(9) is 

/oo 

f(9)e-' iek d9 = f(k) (2.14) 
-00 

with corresponding inverse 

^ [/(*)] = 7T / f( k ) et6h dk = f(°) ■ ( 2 - 15 ) 

J_ 00 



Inserting ( 2.13| ) into the F.T. of ( |2.11| ) (first taking —^e 9 to the left hand side to ensure 



the existence of the transform), we obtain 

Hfdo) + \e e ] = 2, tt V s(k) (1 + Y^y) - 2 • L R {k) (${k) + ) • ( 2 - 16 ) 
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This can be compared with the Fourier transform of the kink limit of fl2.1| ), which is 



2in 2 a5(k) - 2 i L R (k) <p(k) 



(2.17) 



The two will match if 



and 



m = m + 



m 2 

i - 4>(k) 



a 



a 



1 + 



X(0) 



(2.18) 



(2.19) 



1-0(0), 

After some elementary manipulations we can invert the Fourier transform of x, resulting 



m 



X(0) 



3 ikd sinh(^)cosh(^) 



dk 



cosh(^)sinh(^(e-l))2vr 



and, from ( p. 19 ) 



a 



a. 



(e-i) 



(2.20) 



(2.21) 



Since ( 2.18 ) involves only x(k) 2 we have made a choice of sign when taking the 
square root. We took the negative sign since only then do the values taken by c e fr(oo) 
when flowing from a minimal model M p , q always assume the form 1 — Q/p'q', as has to 
be the case if the infrared destination of the flow is to be another minimal model. It 
is possible that the flows found for the other sign choice also have an interpretation in 
some other context, but we will not explore this question here. 



3 The flows 

By construction, the new massless system exhibits the same ultraviolet central charge 
as the corresponding massive system: 

Furthermore, the fact that the kink limit of the new system can be mapped exactly onto 
that of the massive system makes it natural to suppose that the recipe for choosing £ 



and a given by (|2.4D , ( j2.5|) and ( |2.6| ) also holds for the massless perturbations. We take 
this as a working hypothesis for now; in the next section it will be subject to some 
more detailed checks. However, one immediate difficulty should be mentioned: with our 
convention that the coprime pair (p, q) specifying M. Pt q has p < q, the value assigned to 
£ by (|2.4|) for <f>\2 perturbations is less than 1. The same is true for ^15 perturbations 



of M Pt q whenever q > Ap. This is not a problem for the massive perturbations, but 
the massless kernel 4>(9) has a pole at 9 = |(£ — l)7ri, which crosses the real 9 axis as 
£ dips below 1. In addition, the formula ( p. 21 ) clearly has a pole at £=1. All of this 
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leads to technical complications, and, while they could perhaps be overcome by analytic 
continuation, we will leave the issue to one side for now. In the rest of this paper we 
will focus solely on the 02i and 0i5 perturbations, found via equations Q2.5| ) and (2.6) 
respectively, and restrict the 0is perturbations to models M. VA with q < Ap. 

In the infrared, the massless system mimics a massive system with £ replaced by 
£— 1 and a replaced by a' = a£/(£— 1), and so 

CeffM = 1 - ^V) 2 = 1 " 77^. (3.2) 
Consider first the perturbation of Ai pq by 02i- For it to be relevant, we must impose 



that p > q/2. Substituting the requisite values of £ and a (taken from (pjjj)) into (3.1) 
and (p.2|) we find: 



6 6 

M p , q + <p 2 i ■ P>q/2, c efT (0) = l , c cfr (oo) = l-- -. (3.3) 

pq [q-p)q 

Similarly, for the 0is perturbations we have 

6 6 

M p ,g + 015 : P < q/2 , c eff (0) = 1 , c cfr (oo) = 1 - — . (3.4) 

pq p\w—q) 

As promised at the end of the last section, the infrared limiting values of the effec- 
tive central charges are all consistent with the destinations of the flows being minimal 
models^]. Unfortunately, a knowledge of c e fj alone is generally not enough to identify a 
nonunitary minimal model uniquely, so the destinations cannot be completely pinned 
down by this information. Moreover, there are cases in which even the knowledge of 
the full c e ff (r) is not sufficient to identify the model or models involved. The possibility 



of such an ambiguity was first noted, for massive perturbations, in [37|. It goes by the 
name of the 'type II conjecture', and explicitly it equates the effective central charges 
of the following models: 

M Pt g + 0i5 <-> M q / 2 ,2p + 021 <-> - A/ / 2 p,< ? /2 + 012 (p = 2n+l, q = 2m) . (3.5) 
The second equivalence is just a question of labelling, but the first is much less trivial. 



However, all follow from the massive NLIE (2.1) and the recipe ( |2.4| ) ( p76| ) , assuming 
the correctness of these equations (see |25| ). Clearly the massless theories, whose NLIEs 
are parametrised according to essentially the same recipe, will suffer similar ambiguities; 
some examples will be encountered later. 



In spite of these provisos, the forms taken by the results (3.3) and (|3.4j) make for 



some obvious conjectures, and these will receive further support, both analytical and 
numerical, later in the paper. Such checks will never resolve type II ambiguities, so for 
these cases the conjectures are better motivated by the belief that the pattern observed 



%ad the opposite sign choice for x($) been taken in (2.20), the values of c e ff(oo) for the $21 and 



perturbations would have been 1 — 6(q—p)/(qp ) and 1 — 6(4p— q)/(pq ) respectively 
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in other cases should hold in complete generality. Taking into account the l p < q 1 
labeling convention for a minimal model M P , q , the 0is case splits into two and we have 
the following predictions for massless perturbations of general minimal models: 

M p ,g + (f>2l -» M q -p,q ( p < q < 2p) , (3.6) 
M p . q + </>i5 -» M pAp -q (2p<q< 3p) , (3.7) 
M p . q + 015 -> M Ap - q , p {3p<q< 4p) . (3.8) 

There is only a single flow on this list from any given minimal model: when 02i is 
relevant, 0i5 is irrelevant, and vice versa. Note also that the case (|3.8| ) can only occur 
as the last member of a sequence: for 3p < q < 4p, £ lies between 1 and 2 and the next 
step would therefore involve a value of £ less than 1, and we have already decided to 
exclude such cases. 

The previously-known flows (|Fl]) and ( |1.2[ ) are reproduced on setting q = 2p — 1 in 
(|3.6[) and q = 2p + 1 in ( |3.7D . In order to understand the more general pattern, it is 
convenient to define an an 'index' I = 2p — q, and to rephrase (|3.6|) and (|3.7|) as 



AW/ + 02i - AVj,2p-j, (^a') = (f-l ) ^l) , (3-9) 
M P ,2 P+ I + 015 -» A4p,2p-J, (£,«') = (^,2^7) • (3-10) 

For reference we have included the values of £ and a' that should be used in the NLIE 
in order to dial up the corresponding flow. In all cases -f(IR) = — /(UV); the flows (1.1) 



and (|l.2j ) make up the unique sequence with ]/] = 1. For \I\ > 1, there may be more 
than one sequence, with flows sharing the same value of \I\ interlacing each other. This 



is reminiscent of the generalised staircase models studied in 34], but is a little more 
complicated owing to the constraint that the pair of integers labelling a minimal model 
must always be coprime. The number of different sequences with index ±1 is therefore 
given by the Euler 93-function y>(|J|), equal to the number of integers less than |/| which 
are coprime to |/| (so for n = 1 . . . 6, ip(n) = 1,1,2,2,4,2). The index measures the 
distance from the line q = 2p across which the relevance and irrelevance of the fields 
02i and 0i5 swap over. All of this is perhaps best seen pictorially, and in figure |l| some 
of the predicted flows are plotted, superimposed on a grid of the minimal models M p ,q- 
The horizontal arrows of length |J| correspond to 02 1 perturbations, while the vertical 
arrows, of length 2|/|, are 0i5 perturbations. 

As a first check on our results, we evaluated the effective central charge numerically 
for a number of members of the \I\ = 1 series, and made a comparison with the results 



from the massless TBA equations discussed in [19]. In 1£] such equations were written 
in a 'universal' form of the kind first described in [35], but for numerical work it is more 
convenient to write the equations in the following way: 

" />oo 

Sa(0) = Va{e) ~ ~ S *b) / " 6 ') ln (! + e ~ £b{9 ' ] ) d6> » (3- 11 ) 



6=1 
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10 

p 



Figure 1: The grid of minimal models M p , q and a selection of the predicted flows, showing the 
unique sequences with \I\ = 1 and \I\ = 2, and one of the two sequences with \I\ = 3. Also shown are 
the lines q=p, q=2p and q=4p. 



Ceff(r) 



o ^ f'OG 

Va(e)ln(l + e- £ ^) d6 (3.12) 



where n = £ — 3 is integer, is the incidence matrix of the A n Dynkin diagram, 

Va(0) = ^{e e 5 aA + e- e 5 a>n ) , (3.13) 

and 

The agreement between the NLIE and the TBA is extremely good, and is illustrated in 
table [j] and figure [2]. 

Satisfied that our NLIE correctly matches the previously known flows, we can turn 
to the new families with index |/| > 1. Figure ^ shows the (unique) series with |/| = 2, 
and figure || one of the two possible series with |/| = 3. For all of these flows the 
effective central charge initially increases from its UV value, oscillates as the system 
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Model 




a') 


r 


TBA 


NLIE 


M 3 ,7 + 015 


(6, 


!) 


0.01 
0.02 


0.71421869185983 
0.71408433022934 


0.71421869185981 
0.71408433022935 


A4 3j5 + <p21 


(5, 


1) 


0.01 
0.02 


0.59996121217228 
0.59986582853966 


0.59996121217225 
0.59986582853967 



Table 1: Comparison of c e ff calculated using the TBA equations and the NLIE. 



0.8- 




-1 0,1 
log (MR/2; 

Figure 2: TBA (boxes) and NLIE data (solid lines) for the 'diagonal' sequence with \I\ = 1. 



size is increased, before finally flowing to the predicted IR fixed point. This perhaps- 
surprising behaviour is in contrast to the |J| = 1 series where all of the flows are 
monotonia We will return to this point at the end of section 4. 

Note that no TBA equations are known for the massless flows with |7| > 1. This 
puzzle can be highlighted by looking at some cases where the corresponding massive 
TBA system can be conjectured. 

We consider three families of massive systems found in the 'ADET' class of models 
classified in The first set is obtained from Q3.ll ) simply by replacing the driving 
term (gig) by 

v a (9) =r5 aA cosh6 . (3.15) 
This gives the massive flows corresponding to the |/| = 1 massless flows. For the 



second set, we continue to use the driving term (|3.15|) , but replace l^"^ with l^ b n ' (the 
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incidence matrix of the 'tadpole' graph T n = A2 n /^2), letting 1 label the node furthest 
from the 'tadpole' node (see figure ||). In [36| these systems were identified with the 



models M n +2,2n+2 + <t>2\ for n odd and M n +i,2n+A + 015 for n even. These are precisely 
the models in the |/| = 2 series. The massless NLIE predicts, in addition to these 
previously- known massive flows, the existence of interpolating flows (n — * n — 1) within 
this family. 

The final family of massive TBA equations is obtained by replacing the A n incidence 
matrix with the D n one. For these cases the type II conjecture mentioned above plays 
a role, and for each TBA system there are two possible identifications: 

A+ = M 2 n-l,2+4n + 015 and A~ = Mi+2n,4n-2 + 021 • (3-16) 

The sets A^ correspond, in the ultraviolet, to the two series of models with index 
|/| =4. Turning to the massless flows implied by the NLIE, there is a further ambiguity 
in the infrared destinations. However, the general pattern of (|3.6|) , ( |3.7|) and (|4.6|) — ()4.8[) 

suggests 

At^Al_, . (3.17) 

We have checked the low-lying members of each family of massive TBA equations 
against the massive NLIE ([2.1]), finding agreement to our numerical accuracy (about 14 
digits) in each case. However, returning to the question of finding massless versions of 
these equations, we observe a key difference between the first family and the other two: 
as is clear from figure [5l the graphs for the |/| = 1 systems have a Z2 symmetry which is 
in general absent from those for |/| = 2 or 4 (the single exception occurs when |/| = 4, 
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-4-2 2 

log(r/2) 

Figure 4: One of the two sequences with |7| =3. 



n = 4). The trick that allowed us to move between massive and massless systems for 
\I\ = 1 by swapping the z/'s defined in ( 3.15] ) for those defined in ( p. 13 ) relied completely 
on this symmetry. This 'Z2 trick' was first employed in ji; more generally, the nodes 
T and t n' could be replaced by any pair of nodes related by a Z2 graph symmetry. 
To the best of our knowledge, all massless TBA systems that have been discovered to 
date are related to massive systems in essentially this way. Thus if the T n and D n 
TBA systems do have associated massless versions, they are likely to be of a somewhat 
different nature to all previously-encountered examples. 

The \I\ = 4, n = 4 system might appear to offer a counterexample, since the 
exceptional symmetry of its graph does allow a massless version of the massive TBA to 
be constructed. However this equation turns out not to reproduce the massless NLIE. 
For example, for the massless TBA c e ff(oo) = 5/7, which identifies the IR destination 
of the flow as rather than the A^io or found by the NLIE. Furthermore, 

the short-distance expansion of c e s(r) turns out to be inconsistent with the massless 
TBA flow being related to the massive one by an analytic continuation of the coupling 
A - rather, it is the flow produced by the massless NLIE which has this property. This 
failure of the Z2 trick to produce the analytically-continued massless flow can be put 
into a more general context. Recall that TBA equations have associated Y-systems, 
and that these entail a periodicity Y a {9) = Y a {6 + iP) for certain functions Y a (9), 
where P depends on the particular Y-system (see |35| and also jl3]). Suppose that a 
diagram symmetry relates nodes a and a for this system. Then it can be argued that the 
associated massive and massless TBA equations will be related by analytic continuation 
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A n : 



n-1 



1 n : • O O O O 00 



n-1 n 



p n 



D n : 



Figure 5: The A n , T n and D n graphs, indicating in each case the node to which the driving term 
rcoshf in the massive TBA system should be attached. 



if Y a (6) = Ya(9 + iP/2), and not otherwise. In particular, for systems related to the 
ADET diagrams, such a property holds if and only if the nodes a and a are 'conjugate', 
where conjugation acts on TBA diagrams in the same way as charge conjugation acts on 
the particles in an affine Toda field theory |3^.|3^]. For the D n diagrams, the fork nodes 
are related by charge conjugation for n odd, but not for n even. This thus matches 
the observation above that the massive and massless Z^-related systems are not related 
by continuation in A; it also helps to understand from the TBA point of view why the 
models of pi are related to H$ by analytic continuation for N odd, but not for 
N even. 



We conclude this section with some further observations. Following [36| we note 
that two models in the |/| = 2 set possess N=l supersymmetry in the ultraviolet. In 
the notation of 36] we have 



SM 2 , 8 + (f>\f * M 3)8 + 0i5 , (3.18) 

and 

SM 3J + ~ M 7A2 + 02i • (3.19) 

The operators <fy^ and <$° p are SUSY-preserving, so the massless flows M.-&$ — > M.^^ 
and A^7,i2 — > M.5,12 should exhibit spontaneous breaking of N=l SUSY. The situation 
is reminiscent of the flow between the tricritical Ising and Ising models discussed in j|, 
fMj. In particular, the flow ^3,8 — > M.3,4 is to a theory of a single massless free 
Majorana fermion, just as was the case in fi|,|4(|. This particle, in turn, was identified 
in |4,40 1 as the massless Goldstone fermion. Finally, we mention that that the massive 



TBA system for A^3,8 + 4>i5 was alternatively obtained in [36| as a folding of the N=2 
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supersymmetric model with spontaneously-broken symmetry of [41|, for k=3. This 
suggests that there should also be a (probably nonunitary) flow from the 'parent' N=2 
supersymmetric theory into a massless free fermionic theory. Given that a similar 
interpolating phenomenon is also present for ^=20 it is natural to suppose that the 
same behaviour will be exhibited for all of the Z^-related models of [41]. It would be 
interesting to check this conjecture using the aj^j-generalisations of our NLIE. 



4 Perturbation theory 

The claims of the last section can be put on a more secure footing by taking a closer 
look at the behaviour of the effective central charges at small and large r. For this we 
will need the leading asymptotics of the kernels (f)(9) and as 6 — > — oo. That of <f> 
can be found from the residues of the poles in the integrand of fl2.10|) at k = — i and 
k = — 3£/(£— 1), and is: 

^ ™n(f(e-l)) 6 + 7r(£-l)cos(^) + ---' oo.^.i) 

Similarly, 



X(0) 



? e _ 3 S m((^ TT )cos(^ TT ) 

2^sin(f(e-l)) e n^-l)co S (^) + "" °° ' ^ 

We first analyse the NLIE as r —* oo, to see whether its behaviour is compatible 
with the claimed infrared destinations of the massless flows. In spite of the fact that 
conformal perturbation theory about the infrared fixed point is not renormalisable, there 
are a number of unambiguous predictions against which the equation can be checked. 
The key ideas are set out in [§|,|j|], and are further discussed in, for example, || 16,17]. 

In general, the infrared model will be described by an action of the form 



S = + ii\ J if) d 2 x + fi2 J TT d 2 x + (further terms) 



(4.3) 



where is the action, and if) one of the (irrelevant) primary fields, of the infrared 
conformal field theory, ifj might be absent, in which case the only operators attracting 
the flow to the IR fixed point would be the descendents of the identity, TT being the 
least irrelevant example. On dimensional grounds, the couplings /xi and fi2 are related 
to the single crossover scale M as 

Ml = Kl M 2 - 2h , /x 2 = k 2 M- 2 , (4.4) 

with Kx and k 2 dimensionless constants and h the conformal dimension of ifi. 
'for more on the k—2 case, and its connection with the theory of dense polymers, see [po|,^2| 
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Standard methods of perturbed conformal field theory can now be used to calculate 
the first corrections to c e fj(oo), with results that are good at least up to the order at 
which the 'further terms' cut in. The first perturbing term in ( |4.3| ) yields a series in 
fiiR 2 ~ 2h = K\r 2 ~ 2h if all powers of n\ contribute, or fi 2 R 4 ~ ih = k^t 4-4 ' 1 if only even 
powers appear, as happens when tp is odd under some symmetry of 5j* R . The second 
term always contributes a series in H2R~ 2 = K-2 r ~ 2 , the first two terms of which were 
found explicitly in p3j. Putting everything together gives c e fj(r) the following large-r 



expansion: 

c eff( r ) ~ c fj(oo) + (a series in K^r 2 ~ 2 ^ or k 2 7^—4/1^ 

c cs(°°) k 2 r~ 2 H Cc fii®°) K 2 r -4 _|_ ^f^^gj- terms) . (4.5) 

In contrast to the UV situation to be discussed shortly, this series is only expected to be 
asymptotic, and furthermore there is very little control over the omitted 'further terms'. 
Nevertheless, it allows for some useful comparisons with results from the NLIE. 

When r is large, the nontrivial behaviours of the functions Jl and fn appearing 



in (2.7) and ( |2.8| ) are concentrated in 'kink' regions near 9 = — lnr and 6 = lnr 
respectively. So long as the contours C\ and C2 are kept a finite distance from the real 9- 
axis, the functions ln(l+e ± ^ i ) and In^l+e^^) appearing in these equations are doubly- 
exponentially suppressed in the central zone between the kinks, and so the principal 
interaction between equations ( pT7| ) and (pTg|), and hence the principal correction to 
c e g (r), comes from the exponential tail of the kernel function x($); given by ([4.2;). Since 
we are only worrying about the first few terms in the expansion, and we are free to vary 



£ so as to avoid any 'resonance' effects, we can discuss effects of the two terms in (4.2) 
separately. Consider first the second term, decaying as e 3 ^" 1 ) as 9 -> -00. Inserted 
into ( f2.7j ) and considered iteratively, corrections to /r(9) as a series in r~ 6 /(£ -1 ) will be 
generated. Feeding through into c e fr(r), these corrections can be matched against the 
K\ terms in (ff.5|), allowing h, the conformal dimension of the field ip, to be extracted. 
Comparing with the Kac formula h a b = ({bp'—aq') 2 — (p'—q') 2 )/(Ap'q') at the appropriate 
(IR) values of p' and q' leads us to conjecture the following pattern of arriving operators 

M Py q + 4>2i -> M q -p, q ( p < q < 2p) arriving via foi , (4.6) 
M PA + 15 -> M pAp -q (2p<q < 3p) arriving via ^15 , (4.7) 
Mp, q + 0i5 -» M4 P -q,p (3p < q < 4p) arriving via (f> 51 . (4.8) 

In making these identifications, we took account of the fact that 4>2i is an odd operator, 
while is even. 



Exactly the same line of argument starting from the e term in (4.2) reveals a further 
series of corrections to c e g as powers of r~ 2 , perfectly adapted to match the TT terms 



in (4.5). However this time it is possible to say more, by observing that the effect of 



this part of x(0) on equations ( [2.7j ) and ( |2.8[ ) can be reabsorbed into a shift of r by a 
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constant, and furthermore that this constant can be expressed in terms of c(r). This 
allows the first few terms of the iterative expansion to be found exactly. The same idea 
was used in [17| to extract IR asymptotics from various massless TBA systems; in the 



current context we find 



7TC eff (oo) 2 _ 2 27T 2 C eff (oo) 3 _ 4 

sinf(£-l) r sin 2 f(£-l) r 



Ceff(r) ~ Cefr(oo) - ^t— + ^ 2 n '— r 4 + (4.9) 



Matching r terms in ( [4.5D and ( |4.9| ) gives the exact relation K2 = 6/(7r sin^(£— 1)) ; 
the fact that the r -4 terms then agree provides a nontrivial check on the IR behaviour 
of the massless NLIE. 

We also performed some numerical fits on the IR data. Our results are summarised 
in table ||. Accuracy was fairly low, but where relevant we have also included the 
predicted ('exact') values of coefficients, obtained from equation 



Model 




Ceff(r)-c eff (oo) 


M3,5 + <5 
Exact 


!»21 


0.58041577r" 2 + 0.000061 r" 3 + 1.66r~ 4 + 0.12r~l + . . . 
0.58041579 r~ 2 + 1.6844 r" 4 


M 3 j + <f 

Exact 




1.30591 r- 2 + 1.65488 r~f + ... 
1.30593 r- 2 


M.4,7 + (f 
Exact 


!»21 


0.51029 ln(r) r~ 2 - 0.68683 r~ 2 + ... 
0.51020 ln(r)r- 2 


Exact 


021 


-6.8516r~ 2 + 8.8213r"t + . . . 
-6.8510 r- 2 


M 5 fi + <? 
Exact 


!>21 


-1.1253 ln(r) r~ 2 + 1.3757r~ 2 + . . . 
-1.125 ln(r)r~ 2 


M 6 ,n + 
Exact 


015 


-2.32081 r- 2 + 12.85 r' 4 - 1.527 + . . . 
-2.32082 r~ 2 + 12.56 r" 4 



Table 2: Comparison of the infrared expansion of c e fi(r')— c e ff(oo) for various models against the 
exact coefficients. 

The situation in the ultraviolet is in many respects much simpler. Conformal per- 
turbation theory gives direct access to a function c per t(r), related to the ground state 
energy ( |1.3| ) as Eq(M,R) = — ^>c pert (r). Note that c pcrt contains a bulk part which 
must be subtracted before comparisons are made with the NLIE. For a theory per- 
turbed by a relevant primary operator <fi with scaling dimensions (h vv , h vv ) , c pe rt has 
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the expansion [44,45] 



^pcrt 



(r) = c(0) + ^2C n (\Ry) 
n=l 



(4.10) 



and, in contrast to the situation in the IR, the series is expected to have a finite radius 
of convergence. Here y = 2 — 2h uv , A is the coupling, and the coefficients C n are given 
in terms of the connected correlation functions of the perturbing field on the plane as 



12 (- 1) 
nl (2tt) 



yn—l 



d 2 Z 



II|7¥ (^(0)0(1, 1)^2, ^2) 



i=2 



(z„,z„)y(oo)) ( 



(4.11) 



where V creates the CFT ground state on the cylinder. This is the state with lowest 
conformal dimension, and for a general minimal model M Pt q it corresponds to the field 
00 = 4>ab with a and b integers satisfying bp — aq = 1. (Only for the unitary models 
Ai PtP+ i does it coincide with the conformal vacuum (fin.) 

For the perturbing operator, we will be interested in <fi = 1^15 and (fi = 021- As 
mentioned above, ^21 is odd, so only the even coefficients Cm are nonzero for this case. 
Since A must be related to the crossover scale M as 



A = kM v 



(4.12) 



with k a dimensionless constant, ( 4.10| ) is a series in r y for the </>is perturbations, and 
r 2y for the (p2i perturbations. 

The effective central charge calculated using the NLIE is expected to expand as 



Ccff 



(r) =c cS {0)+B(r)+J2 c ^ 



yn 



(4.13) 



n=l 



for some value of y. Periodicity arguments suggest that this will be a series in r 6 ^ 1+ ^\ 
and if £ is chosen according to (2Jj) or (|2.6|) then y = 2 — 2h vv and the perturbative 
expansion ( |4.10 ) is matched, with h vv the conformal dimension of either 1^21 or <fii^, and 
all odd terms zero for the (p2i case [^J. The irregular bulk term B(r) must be subtracted 
before c pert can be compared with c e ff. Fortunately, it can be obtained exactly from 



the NLIE, using a small generalisation of arguments used in [44| and [21]. The term we 
need is given by the behaviour as r — > of 



3ir 
2^2 



c 2 



,d_ 

de 



l Q n +e fL(o)) d9 _ f e -e± ln , 1+e -f L (eu de 

jc! de 

>0 



(4.14) 



where the '> 0' indicates that only those parts of the contours C\ and C2 with positive 
real part should be taken, and the symmetry between /l and fn was used to trade the 
first two integrals in (^^) for the prefactor 2. Consider the r — > limit of ( |2.7| ) and 
(2.8) in the region <C ?R.e0 ^ ln(l/r) , where the driving term —i^e~ e in ( |2,8| ) can 
be dropped. Take the derivative of these equations with respect to 9, and extract the 
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contributions to the convolutions proportional to e using (4.1) and ( |4.2| ). These should 
cancel either against the remaining driving term or between themselves to ensure that 
the functions /t and Jr have no such dependency. This leads to the equations 



r v^sintef) 
i- 
2 



vrsin(f(e- 


1)) 


3 




27rsin(§(£- 


-1)) 


V3sin(f 





7rsin(§(£- 


-1)) 


3 




2vrsin(f(e- 


-1)) 



/ e ^'4ln(l+e^( e '))^'- / e -»'Ab(l + e-/«(»V 



c 2 



9 i d 

do 

cji d 

de 



7 ln(l+e-^ e '))d#' 
7 ln(l+e f L ^)d0' - 



c 2 

X) 



de 



7 ln(l+e /i(fl,) )de' 
7 ln(l+e- /i ( e '))^' 



jc 2 yci do' 



Solving for the integrals needed for the evaluation of ( 4.14 ), we find 

B-massless (0 



3 sm(fQsin(f(g-l)) ^ 

7T 



sin(-7r^) 

A similar if slightly simpler analysis of the massive equation (|2.1[) predicts 



(4.15) 



B 



massive 



(r) 



V3 sin(f 



(4.16) 



2 7rsin(f(C+l)) 

It can be checked that this formula combines the results given separately for ^12, 4>2i 
and (pis perturbations in [46,47]. 

The formula ( 4.15 ) has a pole when £ + 1 is an integer multiple of 3. Since the 
overall result remains finite this infinity should cancel against one of the terms in the 
regular expansion of c pe rt, leaving a logarithmic contribution to the final result (see for 
example [17, 48|, 49|). At this point the calculation might appear to split into two cases: 
the term proportional to r 2 in the perturbative expansion occurs at an order re which 
depends both on the dimension and on the 'parity' of the perturbing operator: n = 2/y 
for , and 1/y for 021 • However y = 6/(£+l) for ^15 and 3/(£+l) for 021 p5| , so 
re = (£+l)/3 in both cases and they can be treated simultaneously. The logarithm is 
found by evaluating 



lim 3sin(|Qsin(|(e-l)) , 2 



£^3n-l 



7T sin(7r£) 



2n 



M-l 



which yields 



B 



massless 



£=3n-l 



2vr 2 n 



r 2 In r . 



(re G N) 



(4.17) 



(4.18) 
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Similarly, and for exactly the same values of £, logarithms arise in the massive pertur- 
bations. In these cases we find 



B r , 



(r) 



£=3n-l 



2vr 2 n 



r 2 lnr 



In G 



(4.19) 



This seems to leave the coefficient of r 2 term unknown, but one piece of information 
can be extracted, which will be used later. Recall that the logarithm is linked to 
a divergence in the perturbative integral defining the term C n . Providing the same 
regularisation scheme is used for both flow directions, the perturbative contributions 
should cancel upon taking their sum or difference. The remaining finite contribution 
can be found directly from ( [4.15| ) and ( 4.16 ): 



Bmasslessy) i Bmassivey) 



£=3n-l 4-7T 



(n e 



(4.20) 



where the plus signs should be chosen if n is odd and the minus signs if n is even. 



Using an iterative method to solve (2.7) and ( |2.8| ), c e s(r) can be calculated to high 
accuracy. After subtracting the ultraviolet central charge and the relevant bulk term, 
the perturbative coefficients c n can be estimated via a polynomial fit. To make the 
comparison with CPT, we must also fix the value of k. For the massive systems this 
was determined exactly in ^q, ^7j, and it turns out, numerically at least, that the 
same relationship between A and M holds in the massless cases, apart from a factor of 



either — 1 or i. As explained in [19], the relationship between the massive and massless 
perturbations depends on the parity of the perturbing operator. For 0is , the massive 
behaviour is related to the massless by flipping the sign of the coupling constant A. This 
has no effect for the odd operator 021 ; instead the required transformation is A — > iX. 
The CPT coefficients C n are the same for both flow directions so, provided the mass 
and the crossover scales M are equal, we expect 



»15 

hi 



C-n — ( 1) C n , 
C2n = (-1)" C 2n , 



(4.21) 

(4.22) 



where c n denotes the expansion coefficients obtained using the massive NLIE. Tables ^ 
and U present the first few perturbative coefficients for the models M^u and Ai^s 
perturbed in both massless and massive directions. The relative signs of the coefficients 
confirm relations (4.21) and (4.22). Similar results were obtained for the models , 
-^5,12 and A^3,io perturbed by ^15 , and the models , ^3,5 and ^7,11 perturbed by 
(p2i- These provide strong support for our claim that, modulo the bulk terms, the results 
from the massive and massless integral equations are related by analytic continuation. 

For the (^15 perturbations, it is also simple to check the value of k directly, since the 
first term in the perturbative expansion Q4.10 ) is just given by 



d = -12(27r)( 1 -^C, 



00 015 00 



(4.23) 
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n 


c n (massive) 


c n (massless) 


Acn 







19/22 


19/22 


1 




1 


-0.393148695201 


0.393148695201 


±7- 10- 


-13 


2 


0.0166115829 


0.0166115834 


±5 • 10" 


-10 


3 


-0.001225836 


0.001225827 


±8 • 10- 


-9 


4 


-0.00003425 


-0.00003422 


±2 • 10- 


-8 


5 


0.0000116 


-0.000012 


±3 • 10- 


-7 



Table 3: Comparison of massive and massless UV coefficients for A^4,ii+0ib- The last column 
reports the difference between the absolute values of the two. 



n 


c n (massive) 


c n (massless) 


Ac n 







0.85 


0.85 


1 




2 


-0.1128538069088 


0.1128538069085 


±3 • 10" 


-13 


4 


0.1607667041 


0.1607667048 


±7-10" 


-10 


6 


0.003799534 


-0.003799535 


±1 • 10" 


-9 


8 


-0.0009586 


-0.0009585 


±1 • 10" 


-7 


10 


-0.0001554 


0.0001554 


±8 • 10" 


-8 



Table 4: Comparison of massive and massless UV coefficients for A^5,8 + <^2i- The last column reports 
the difference between the absolute values of the two. 
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Model 


num 


K 


M 5 ,12 


0.012665147953 


0.012665147955 . . . 




0.0247653386712 


0.0247653469711 . . . 




-0.06407485530 


-0.06407485531 . . . 


M7,n 


0.0092588732986 


0.0092588732985 . . . 


■M. 10,23 


0.02340788397 


0.0234078837 . . . 


M 3 ,w 


-0.0188367 


-0.0188368 . . . 



Table 5: Comparison of exact and numerical values of k 2 for ^15 perturbed models. 



where C, 



00 015 00 



, the operator product coefficient between the perturbing operator ^15 



and the ground state 0o , can be found in |30]]. Comparing ( 4. lCj ) and ( 4.1 3|) at order 
n = 1, we have 

* = ^ . (4.24) 



Using this formula and the massless NLIE, we estimated k 2 numerically for a number 
of models. Table || reports the results, and compares them with the exact expression 
for n 2 for (massive) (^15 perturbations given in [47|: 



4 2 (£ - l) 2 7 ( 



2+e \ -.( sg 



vr 2 (2-3C) 2 (2-0 2 7 2 (fr|) 



v^r(i) r(|) 



12 



(4.25) 



where 7(2;) = r(i-i) anc ^ ^ ^ s ^ ne mass °f the lightest kink.[] The agreement is clearly 
very good. Note that k 2 is sometimes negative - these are cases where, in the normali- 
sations of 5C ] , Cfo 16 turns out to be purely imaginary. 

Now we return to the behaviour of the flows (|lJ) and ( |1.2[) , to mention one reason 
why, exceptionally, c e ff(r) should be a monotonic function of r for these flows. The 
behaviour of a flow is determined at small r by the first nonzero term in the perturbative 
expansion. Thus whether c e ff initially increases or decreases will typically be determined 
by the sign of C\\ (</>is) or C2A 2 (^21)- In particular, if C\ (respectively C2) is nonzero 
then for one sign of A (or A 2 ), c e ft will initially increase, leading to an immediate violation 
of the 'c e ff-theorem', and a flow which must be non-monotonic if c e fj(oo) is to be less 
than c e ff (0). But for the model -A^s + 4>2i, the first coefficient C2 was calculated to be 
zero in Jl9[ . In this case the asymptotic UV behaviour is instead controlled by C4A 4 
and this permits c e ff to decrease initially for both (massless and massive) signs of A 2 . A 
similar calculation for the other models in the |/| = 1 series finds C\ ((^15) or C2 ((^21) 
to be zero and thus, as in the first model, all of the flows are able to be monotonic. A 
numerical fit of the data for models higher up in the series confirms these coefficients are 
zero within our numerical accuracy. For all other sequences of models, such vanishings 

"note the relation was given in @ in terms of £ FLZZ = ^ and m, the mass of the lightest breather, 
related to M as m — 2M sin(f£) 
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of C\ or C2 do not occur, forcing at least one of each pair of massive and massless flows 
to be non-monotonic. One remaining mystery, out of many, is why it should always be 
the massless flow which is non-monotonic. 

The |/| = 1 sequence of flows was found in [12] via an associated staircase model, 
and the final step of this staircase interpolates in the infrared to a massive model with 
c eff = 0. It turns out that the NLIE also reproduces this behaviour: even though 
the operator 0i5 is not a member of the Kac table for M-2,5 , (|3.7|) formally predicts 
a further flow, to a theory with c e fj = 0: M.2,5 + 015 — ► -M-2,3- There is nothing 
to stop us using the massless NLIE ( ^?7| , |2.8|) and the 0i5 recipe ( |2.6[) to compute an 
effective central charge for this putative flow. Our numerical results show c e g to have an 
exponential behaviour in the infrared, as would be expected for a massive, rather than 
a massless, flow. Furthermore, comparing c c fr (massless) with that of the massive flow 
-^2,5 + 015 calculated using (^) and (^), we find the effective central charges match 
exactly. They also turn out to coincide with the results from the more standard massive 



flow A^25 + 012 7 computed using (2A) and (|2.4j ). This means that the flow is at least 
physically reasonable, since 0i2 is the single relevant primary field in A^s- However, 
it remains a curiosity that the same flow can be found from three different nonlinear 
integral equations. A sample of our numerical results is shown in table ^. Note that 
this example shows that the straightforward prediction of scaling dimensions based on 
periodicity arguments is not always correct: for M.2,5 + 015 one would have expected a 
series in r 6 / 5 for c e fj(r), whereas in fact the expansion is in powers of r 12 / 5 . 



r -^2,5 + 0i5 (massless) M.2,5 + 015 (massive) M.2,5 + 4>i2 (massive) 
~rl01 0.39997512539833 0.39997512539839 0.39997512539863 

0.2 0.39254149935036 0.39254149935037 0.39254149935054 

Table 6: Comparison of effective central charge for the three a priori different models. 



If we compare the infrared destinations of (3.7) and ( |3.8[ ) we see it is possible for 
two different ultraviolet models, both perturbed by 0i5, to flow to the same infrared 
fixed point, one attracted via the irrelevant operator 0i5 and the other via 0si. Figure ^ 
illustrates two such flows. As was noted in §3, the sequence attracted by 4>5i necessarily 
stops at this model, but the other sequence may continue to flow down further. 



For some models, the predictions (4.6)-(|Q|) must be treated with caution, as we 



were unable to check them explicitly. Examples are the minimal models Mp,q perturbed 
by 015 when p and q are related as 4p — q = 1. For these cases the twist a' is equal to 2, 
and a trivial shift of Jl/r(^) leaves a system with a' = 0. This prevents us from tuning 
in a straightforward way onto the required ultraviolet model: numerically, the massless 
NLIE simply recovers a flow from c c fj(0) = 1 to c e fj(oo) = 1 . 

We also observe that ( |3.6D formally predicts the following flows from unitary minimal 
models: 

AW + 021 -» M 1>p , (£, a') = (^r, 1) • (4.26) 
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log(r/2) 
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Figure 6: Two different flows to the same IR point, Mt^ : M$,29 + $15 attracted via </>5i , and 
A^7,ig + 0i5 attracted via 0i 5 . 



However it is a little hard to decide what is meant by -Mi iP , and indeed, while our 
numerical results indicate that these flows behave as expected for small r, at some 
intermediate scale c e g(r) appears to have a discontinuity. This could be linked to a' 
being equal to 1. In any event, we suspect there may be a square root singularity 
though a detailed investigation will have to be left for future work. Note that by the 
type II equivalence, a similar phenomenon occurs for the models M p ^ p -2 perturbed by 
015. In one case a detailed check on these exceptional flows can be made with relative 
ease: the Ising model M^A perturbed by (j>2i- For real coupling A, an exact expression 
for the (massive) effective central charge was given in |]45| . Using this we can make 
an exact prediction for c e g- in the massless case: we send r — > i r in the perturbative 



expansion, swap the logarithmic bulk term for that of the massless model ( 4.15 ) and 
find the coefficient of the r 2 term using ( f4.20| ). The result is: 



3r2 
2tt 2 

6 



lnr 



hl7T + j E + 



7T 



k=l 



2^/3 
' -{2k 



1)tt + 



2{2k - l)vr 



(4.27) 



where = 0.57721556... is the Euler-Mascheroni constant. This formula predicts 
a square root singularity at r=7r, preventing a smooth interpolating flow to the far 
infrared, but we can at least match against the NLIE for values of r out to this point. 
As table f| shows, our numerical results agree well with the exact formula (4.27). 
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0.0001 


0.500000014242146 


0.500000014242144 


0.6 


0.535668551455508 


0.535668551455511 


1.1 


0.49940046725090 


0.49940046725085 


1.4 


0.4131296597211 


0.4131296597212 


1.8 


0.18775453205 


0.18775453226 



Table 7: Comparison of exact and numerical values for c c g for the <f>2i perturbation of A^3,4- 



5 A remark on 13 perturbations 

In this short section we comment on some features of massless 0i3 perturbations which 
are revealed when similar methods are used. In [^] numerical results were only pre- 
sented for the zero twist (c e ff (0) = c e fj(oo) = 1) case, but the same equation can describe 
a minimal model M p . q perturbed by 0i3. One simply has to set the kernel parameter £ 
(p in |^]) to p/(q—p), and the twist a in that paper to ir/p. 

Defining an index J = q — p H, C, = pj J and the flow from £ to £ — 1 is 

M PtP+ j + 0i3 -> Mp-j^p attracted via 3i . (5.1) 

Again, the number of sequences for each J is given by the Euler ^-function, (p(J) ■ As 
(2) 

in the a 2 case, we cannot access all possible minimal models via this equation. The 



kernel (f)(6) in the massless equation of |2Q] has a pole at in(C — 1) which crosses the 
real 6 axis as £ falls below 1. As before, this could probably be overcome using analytic 
continuation but for now we choose £ > 1, requiring 2p > q, to prevent such problems 

occuring. We also find a set of models, A4 P; 2 P -i + 0i3, where a' = 1 and the would-be 

(2) 

IR model is Aii >p . As in the a 2 case, c e fj suffers some sort of discontinuity at an 
intermediate scale for these flows. 



Solving the massless equation in [2C] for the unitary cases (J = 1), we find that 
results from the TBA equations of Q are matched to high accuracy^*! The results for 
J > 1 are perhaps more interesting, since TBA systems for these nonunitary flows are 
not known. Just as for the 021/015 flows with |/| > 1, the monotonicity property is 
lost. Typical cases are illustrated in figures ^ and |8|. 

6 Conclusions 

We hope to have demonstrated the rich structure of flows that can be studied by means 
of a relatively simple nonlinear integral equation. A number of unexpected features 
have emerged, most notably the consistent failure of the effective central charges to be 
monotonic functions of the system size. Indeed, looking at the full set of minimal models, 



**we understand that this has already been checked by Alyosha Zamolodchikov |51 
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log (r/2 ) 

Figure 7: One of the J = 2 nonunitary ^13 perturbed flows calculated using the NLIE of pCj 



both unitary and nonunitary, we see that for massless flows a monotonic behaviour of 
c e g(r) is very much the exception rather than the rule. 

As for future work, the following directions seem natural: 

(2) 

A more detailed study of both the massive and the massless a 2 -related nonlinear 
integral equations is warranted. This should include a more detailed look at their 
analytic properties, and further comparisons with perturbed conformal field theories. 
In addition, the equations should be generalised to describe excited states. 

The type II conjecture remains a curious observed fact, and a deep understanding 
seems still to be lacking. In cases suffering from this ambiguity, a knowledge of c e s(r) 
will never suffice to disentangle the possible destinations of the flows. It would be 
reassuring to confirm that these models behave as conjectured by some other method. 



The work of [20] was based in part on a massless S-matrix. Here we avoided this 
aspect, but it would be interesting to find an S-matrix description of the new flows. 
One might also hope to find integrable lattice models which would yield the massless 
flows in their continuum limits. 

As discussed in §3, we were unable to treat some flows, due to singularities in the kernel 
4>{0) crossing the integration contour. It would be worthwhile to study the NLIE in 
these zones where £ falls below 1, perhaps by analytic continuation. It would also be 
interesting to study the NLIE of pfj] for Q < 1, where a similar difficulty is encountered. 
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of 0. 



One of the normnitary 0i3 perturbed flows with index J — 5 calculated using the NLIE 



The method described in §2 should enable a massless nonlinear integral equation to 
be obtained from any of the (massive) ADE-related systems described in [52,53]. It 
would be interesting to know whether these share the same strange properties that we 



(2) 

have observed here in the examples related to and a 



As mentioned in §4, the initial discovery of the flows and (|1.2| ) in [12] came via a 
staircase model. The new flows also fit naturally into staircase-like patterns, and it is 



natural to ask whether the existing set of known staircase models 12, 34, 54 -5£] could 
be enlarged so as to include these cases. 

• Finally, it remains an open question to find TBA equations describing the non-unitary 
0i3 flows, or any of the new <p2i or </>i5 flows. Why is there no TBA for any non- 
monotonic case? And what has this to do with the monstronP), if anything? 
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